Galactic tide and secular orbital evolution 
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Abstract Equation of motion for the galactic tide is treated for the case of a comet 
situated in the Oort cloud of comets. We take into account that galactic potential and 
mass density depend on a distance from the galactic equator and on a distance from the 
rotational axis of the Galaxy. Secular evolution of orbital elements is presented. New 
terms generated by the Sun's oscillation about the galactic plane are considered. The 
inclusion of the new terms into the equation of motion of the comet leads to orbital 
evolution which may be significantly different from the conventional approach. The 
usage of the secular time derivatives is limited to the cases when orbital period of the 
comet is much less than i) the period of oscillations of the Sun around the galactic 
equator, and, ii) the orbital period of the motion of the Sun around the galactic center. 
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1 Introduction 



Global galactic gravitational field influences motion of a comet in the Oort cloud in 
the form of the galactic tide. The motion of the comet with respect to the Sun is 
important in better understanding of the Oort cloud. This paper presents equations 
for secular evolution of orbital elements for the comet under the gravity of the Sun and 
the galactic tide. We consider equation of motion derived in Klacka (2009a). Results of 
our paper reduce to the results obtained by Klacka and Gajdosi'k (1999) and Fouchard 
et al. (2005) when several physical terms are ignored. The results are compared with 
detailed numerical solution of the equation of motion given by Komar et al. (2009). 
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2 Equation of motion 

We are interested in motion of a comet with respect to the Sun. The comet is in the 
Oort cloud and wc want to describe the cometary evolution in terms of secular evolution 
of comet's orbital elements. 

The Sun is moving at a distance Rq = 8 kpc from the center of the Galaxy. Cur- 
rently, the Sun is situated 30 pc above the galactic equatorial plane {Zq = 30 pc). 
Besides rotational motion with the speed {A — B) x Rq (where A and B are Oort 
constants) the Sun is moving with the speed 7.3 km/s in the direction normal to the 
galactic plane. Positional vector of the comet with respect to the Sun is r = 77, Q. 
Equation of motion is taken in the form 



cf£, _ GMq 



i + {A- B)[A + B + 2Acos{2 ujQt)] £, 



dt'^ r 

- 2A{A ~ B) sin (2 a;of) V 

+ 2{A - Bf (A - r2Zl) Ro Zo cos (ojot) C 



dt2 



dji _ _ GMq _ ^ _ ^.^ ^ 

+ {A - B)[A + B ~ 2Acos (2 wq^)] V 

- 2{A - Bf (A - r2Zl) Ro Zo sin {cjot) C , 



4: n G Q Zo [cos (uJo't) C ^ sin (uJot) v] 



^ = - [4nGo + 2(A' - B')] Zo, 



r= sje + + C2 , 
^0 = ^-5, (1) 

where G is the gravitational constant, Mq is the mass of the Sun and 

A = 14.2 km s"^ kpc"^ , 

B = -12.4 km s"^ kpc~^ , 
A = 0.124 kpc"^ , 
r-2 = 1.586 kpc"^ , 

Q = 0.130 Mq pc"^ , 

Q = -0.037 M0 pc"^ kpc"^ , (2) 

see Eqs. (26)-(27) in Klacka (2009a). If one wants to use other values of the Oort 
constants A and B, he can use Eqs. (22) in Klacka (2009a): 

Q = Qdisk + Qhalo > 

Qa^sk = 0-126 M© pc"^ , 

Qhalo = {'iTTGy''^ [X {Galaxy) + X (disk) + X (bulge)] , 
X{Galaxy) = -{A -B)x {A + 3B) 
X{disk) = - 396.90 km^ s"^ kpc"^ , 

Xibulge) = — 0.65 km^ s~^ kpc~^ , (3) 
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Fig. 1 Two orbital evolutions of a comet with initial semi-major axis 10000 AU in the Oort 
cloud under the influence of gala<;tic tide. Evolutions are obtained from numerical solution of 
system of differential equations given by Eqs. (13)-(17) with the new terms (solid line) and 
without new terms (dashed line). 

The value X{Galaxy) = 611.800 W s"^ kpc~^ holds for A = 14.2 km kpc"^ 
and S = - 12.4 km kpc~^ Eqs. (22) of Klacka (2009a) can be used. 

3 Secular changes of orbital elements 

Perturbation equations of celestial mechanics yield for osculating orbital elements (o - 

somi-major axis; e - eccentricity; i - inclination (of the orbital plane to the reference 
plane - galactic equatorial plane); f2 - longitude of the ascending node; ui - argument 
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Fig. 2 Time evolution of eccentricity when initial eccentricity is close to zero. 



of perihelion; O is the position angle of the particle on the orbit, when measured from 
the ascending node in the direction of the particle's motion, O = u) + f): 



de 

dt 
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6 V ^ 



cos / + 
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1 + e cos / 
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2 + ecos/ 
Fr cos f -Ft — 7 sm / 



1 + e cos / 



_ sin O 

r N — — r cos I 



f V A* P 

p/{l + ecos/) , 

a (1 - e^) , 
GMq , 



„ sin 6> 

-rjV —. — r COS? 



(4) 



where Fr, Ft and Fjv are radial, transversal and normal components of the disturbing 
acceleration. We use — /i e/j / as a central acceleration determining osculating orbital 
elements if we want to take a time average (T is time interval between passages through 
two following pericenters) in an analytical way 
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VIJ_ [ 
a3/2 2nJo 



2tt 2 
9{f) df 



^ ^ 9{f) df , (5) 



a2 vT^T^ 2n 

assuming non-pseudo-circular orbits and the fact that orbital elements exhibit only 

small changes during the time interval T; the second and the third Kepler's laws were 
used: df/dt = -^/MP - conservation of angular momentum, /T^ = /x/(47r^). 
Rewriting Eq. (1) to the form 



d^^ GMq 



d\ _ GMq 



r = + r?2 + ^2 ^ 

Fx = (A- B)[A + B + 2A cos (2 wq*)] ^ - 2A{A-B) sin (2 uot) r] 

+ 2(A - Bf (A - r2Z^) Ro Zo cos {uiot) C , 
Fy = - 2A{A - B) sin (2 cjo*) C + {A- B)[A + B - 2A cos (2 wq^)] V 

- 2{A - Bf (A - r2Z^) Ro Zo sin (ojot) C , 
-Fz = - C - 4 TT G Zo [cos (wo^) C - sin (ujot) ??] , 
Zo = K sin(wz t + ipo) , 
Zo(t = to) = 30 pc = 6.188 X 10^ AU , 
Zo{t = to) = 7.3 km s~^ = 1.540 AU yr"^ , 



y/4-K G e + 2{A^ - 52) 



a;o = ^ - S , (6) 

where to denotes the current time moment and ipo the initial phase, we can find the 
required components Fr, Ft and F/v of the disturbing acceleration: 

FR = F-eR = Fx + Fy cr y + Fz cr ^ , 

Ft = F ■ BT = Fx ex X + Fy bt y + Fz bt z , 

Fn = F ■ Bff = Fx X + Fy bn y + Fz bn Z , 

Br = (cos i7 cos© — sin 12 sin© cosi, 

sini? cos© + cos fi sin© cos i, sin© sini) , 

Bt = (— cosi? sin© — sini? cos© cosi, 

— sini? sin© -|- cosi? cos© cosi, cos© sini) , 

ejv = (sini? sin i, — cos J? sin i, cosi) . (7) 

Inserting Eqs. (6) into Eqs. (7): 

Fr = {{A-B)[A + B + 2Acos{2u>ot)] r (cos 12 cos© — sini? sin© cosi) 
— 2A{A — B) sin {2 uot) r (sini7 cos© + cos J? sin© cosi) 
+ 2{A - B)^ (A - AZo) Ro Zo cos (wot) r sin© sini} x 
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X (cos Q cos O — sin Q sin O cos i) 

+ ^— 2A{A — B)sm(2 u^t) r (cos J7 cos 6* — sin J7 sin© cosi) 

+ {A - B)[A + B -2A cos (2 uot)] r (sin i7 cos + cos i7 sin cos i) 

— 2{A-Bf (A -r2Zo) Ro Zo sin(cJot) r sin6» sini} x 
X (sin i? cos O + cos J? sin cos i) 

+ {— r sin© sini 

— 4 TV G g' Zq r ^coa (oJQt) (cos O cos — sin sin cos t) 

— sin (wot) (sin i? cos + cos f2 sin cos i)] } x sin sin i , (8) 

Ft = {(^ - B) [A + B + 2Acos (2 wo*)] (cos 12 cos — sin sin cos i) 

— 2A{A — B) sin (2 wot) r (sin O cos + cos O sin cos i) 
+ 2{A-Bf (A -r2Zo) Ro Zo cos(woi) r sin0 sini} x 
X (— cos sin0 — sin i7 cos0 cosi) 

+ | — 2^(^ — B) sin (2 Wot) r (cos J7 cos0 — sin sin0 cosi) 

+ (^ - B) + S - 2ylcos(2 Wo*)] r (sini? cos0 + cosi? sin0 cosi) 

— 2{A - Bf (A - r2Zl) Ro Zo sin (wo*) r sin0 sini} x 
X (— sin ^2 sin + cos i? cos cos i) 

+ I — Ct;^ ?" sin sin i 

— 4: n G g' Zq r [cos (wpt) (cos il cos — sin il sin cos i) 

— sin (wot) (sin 12 cos + cos i? sin cos i)] } x cos sin i , (9) 

Av = {{A — B)[A + B + 2A cos (2 wo*)] r (cos J? cos — sin O sin cos i) 

— 2 A{A — B) sin [2 iJot) ''' (sin 12 cos0 + cosi? sin0 cosi) 

+ 2{A~ Bf' (A - A^o) -Ro Zo cos(a;ot) r sin0 sini} x (sin r2 sini) 

+ |— 2A{A — B)sm{2 Wot) r (cos 12 cos0 — sini? sin0 cosi) 

+ {A- B)[A + B-2A cos (2 wot)] r (sin i? cos + cos i? sin cos i) 

— 2(yl - S)^ (A - /2Z^) Ro Zo sin (wo*) r sin0 sini} x (—cosi? sini) 
+ {— ^ sin0 sini 

— 4 -n G Q Zo r [cos [wot) (cos i? cos — sin i? sin cos i) 

— sin (wot) (sini? cos0 + cosi? sin0 cosi)]} x cosi . (10) 

3.1 A method of averaging 

Inserting Eqs. (8)-(10) into Eqs. (4) and making time averaging represented by Eq. (5), 
we need also the relation between the time t and the true anomaly /. The relation dt 
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Fig. 3 Comparison of two orbital evolutions of a comet with initial semi-major axis 10000 
AU obtained by numerical solution of Eqs. (1) (solid line) and Eqs. (13)-(17) (dashed line). In 
both numerical solutions are the new terms included. 



= (r'^/y/Wp) df yields 

if we take f{t = r) = 0. 

Another possibility is to use the Kepler's equation t — T+{a^^'^ /^/JT){E — e sin E) 
[see also Eq. (11) together with tan(£/2) = ^'(l - e)/(l + e) tan(//2)]. Then, instead 
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of Eq. (5), we obtain 



1 

(5> = ^ j g{E) (1 - e cosiS) dE , 



. a Vl - . 

smr = smE , 

r 

cos / = — (cos E — e) , 
r 



r = o (1 — e cosE) , 

t = r + (o^/^/Vm)(.B - e sinB) . (12) 



3.2 Another method of averaging 

The procedure of averaging may not take into account time evolution in uot and ojzt 
during one orbital period if 27r/a;o 3> T and 2n/ijJz ^ T. 

The secular evolution of orbital elements (and other quantitites) can be calculated 
on the basis of Eqs. (4)-(5) and (8)-(10). We define f2 as O = O + wot- The final result 
then can be summarized in the form 

/da\ 2 [p V 7 ■ ■ o 

{ — ) = — a A — Xa Zq sm I cos U , 
\at / Y M 

Xa = 2{A - Bf (A - TaZ^) Ro + A tt G g' , (13) 

/ de\ _ 5ae 
\dt/ ~ ~T 



^J^{4A{A-B) [cosi cos2a; sin2i? 



+ sin 2a) (cos^ O — cos^ i sin^ 17)] 

+ [4 TT G Q + 2{A^ - B"^) - {A - Bf] sin^i sin 2a;} 

+ 5ae J- {A- Bf (A - r2Zo) Rq Zq sini sina; 



X (cosi cos a; sinJ7 + sinai cos/?) 

, 5ae /p" ^ / 7 . . 
H — — . — 4: -K G Q Zq Sim cos w 
2 y /X 



X (cos oi cos i? — cos i sinoi sin i?) , (14) 

= - A{A - B) a [5 (cosi sin 2a; sin 2/2 

— 2 cos 2a; cos^ Q) + (5 cos 2a; — 3) 

X (cos^ Q H 5—^ cos^ i sin^ J?) 

^ 1 — e'' ' 

H — E2£A^ C5g2 gjjj 2a; sin cos Q 
1 — ^ 

— 2cosz sin^ J?)] 



+ 4( 



[4 ^ G , + 2(^2 - B^) - {A-Bf 
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Fig. 4 Comparison of two orbital evolutions of a comet with initial semi-major axis 50000 
AU obtained by numerical solution of Eqs. (1) (solid line) and Eqs. (13)-(17) (dashed line). In 
both numerical solutions are the new terms included. 



X [(sin^ i — e^)(5 cos 2a; — 3) + 2 cos^ i] 



1-e 



~ J-,{A- Bf (A - TaZg) Ro Zq sini sinw 



X [cosi sinw sini?— (1 — e^) cosw cos J?] 
A'kG Q Zo 4— \f> (si 



2(1 -e2) Y 
X (cos w COS Q — cos i sin a; sin ti) 



sm^ i — e^) sinw 
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+ (1 — e^) cosi sin J?] , (15) 



— cosi (2 — e^(5cos2c<j — 3)) sin^ Q] 

- [4 TT G Q + 2{A^ -B'^) - {A-Bf] cost [2 - e^(5cos2w - 3)]} 

+ ^■^^ " ^° ^° 

X (l — + 5e^ sin^ w) sin i7 

+ 2(13^ ^ - 4 ^ G , Zo -- [(1 - e ) sm 12 cos z 

— 5e sinw (cosw cos — cos i sinw sin/2)] , (16) 

f ) ^ sin. {5e^ [4A(A-5) sin^ 

X (cos 2uj cos — cos i sin 2ctj sin O) 

- (4 TT G e + 2(/l^ -B^) - (^- B)^) cosi sin 2a)] 
+ 2yl(yl- ii)(2 + 3e^) sin 2^2} 

1 e y M 
X sin oj cos w sin Q 

~ o\ 4 - 4:n G e Zq cos i [(1 - e^) cos i? 

2(1 -e^) y fi L 

+ Se^cosw (cosw cos i? — cos i sinw sini?)] . (17) 

Eqs. (13)-(17) yield for the a— component for angular momentum per unit mass Hz 
= \J a (1 — e^) cos i: 

(^^"j = - A{A - B) ^{5e^[(l + cos^i) cos2w sin2r2 

+ 2cosi sin 20; cos 212] + (2 + 3e^) sin^ i sin2i?| 
- {A- Bf (A - r2Zl) Ro Zo sini 
X [(1 — e^) cosi cos /2 + 5e^ sin cj 

X (cosi sin w cos 12 + cos a; sin 12)] . (18) 

Eqs. (13)-(17) show secular time derivatives of orbital elements of a comet during 
one orbital period caused by the perturbation acceleration given in Eqs. (8)-(10). The 
secular time derivatives, represented by Eqs. (13)-(18), are generalization of the results 
obtained by Klacka and Gajdosik (1999) and Fouchard ct al. (2005). 

The secular orbital evolution holds if the method of averaging is acceptable. The 
period of revolution T must fulfill the condition (uz + wo) T" <S 2 tt or equivalent 
condition T <^ { 1 / Tz + 1 / To )~^. This is a consequence of the terms Zq sin(a;ot) 
and Zq cos(a)o^)- 
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Fig. 5 Two orbital evolutions of a comet with initial semi-major axis 50000 AU obtained 
from numerical solution of system of differential equations given by Eqs. (13)-(17) with the 
new terms (solid line) and without new terms (dashed line). 



4 Discussion 

Eqs. (13)-(17) produce identical orbital evolution after each of the following transfor- 
mations: 

1. — > + TT, 

2. — ^ — Zq, W — ^ W + TT, J? — > i? + TT, 

Zq ^ - Zq, n ^ n + -K. 

The first transformation represents a symmetry of the Oort cloud under the action of 
the galactic tide. The second transformation holds due to symmetry of the galactic 
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potential with respect to the galactic equatorial plane. The transformation a; — > + 
TT and J? — ^ J? + TT is equivalent to the transformation i — > — i. The third transfor- 
mation can be obtained from (simultaneous) composition of the first and the second 
transformation. 

In order to obtain the opposite signs in time derivatives of the orbital elements, it 
is sufficient to use the following transformation in Eqs. (13)-(17): 
w— >-7r — W, J7— ^TT — i?. 

This transformation represents an antisymmetry of the Oort cloud under the action of 
the galactic tide. 

The total secular time derivative of semi-major axis is not equal to zero. This can 
have a close relation to the result represented by Eq. (31) in Klacka (2009b). Non- 
zero value of secular time derivative of semi-major axis is caused by new terms in our 
equation of motion. The new terms contain Fi , /2 and g' quantities. Namely, the term 
proportional to A — A-^o ™ ^ and r] components of the acceleration and the term 
proportional to g' in component of the acceleration. Two orbital evolutions of a comet 
obtained by numerical solution of Eqs. (13)-(17) are shown in Fig. 1. Both evolutions 
depicted in Fig. 1 have equal initial conditions. Initial values of the orbital elements 
are ai„ — 10 000 AU, Bin — 0.4. oi,,, = 0, — and tin ~ 90°. The Sun is located 
at distance 8 kpc from the galactic center, Zq{0) = 30 pc and Zq{0) = 7.3 km s^^ at 
the time t = 0. Evolution depicted by a black color is calculated using secular time 
derivatives from Eqs. (13)-(17) which have A = 0.124 kpc"^ r2 = 1.586 kpc~^ and 
g' = - 0.037 Mq pc~^ kpc~^. Evolution depicted by a gray color is calculated from 
Eqs. (13)-(17) without the new terms, i.e. by putting A = 0, A = and g' = in Eqs. 
(13)-(17). Semi-major axis of the comet is not constant for the evolution depicted by 
the black color in Fig. 1. Oscillation in semi-major axis depicted by the black color is a 
typical behavior of semi-major axis found from Eqs. (13)-(17) with the inclusion of the 
new terms. The semi-major axis oscillates around the value close to aj„. We did not 
find evolution of semi-major axis with a tendency to a monotonic increase or decrease 
in time. Only oscillations in evolution of the somi-major axis existed. Evolutions of 
other orbital elements are not significantly affected by inclusion of the new terms. 

Eqs. (14) immediately show that e = if ei„ = 0. We found that if ej„ is close to 
zero, then the eccentricity can increase to a value close to 1. If all other initial orbital 
elements are fixed, then the time span needed for the increase of eccentricity is usually 
longer for the comet with smaller value of ei„. Such situation is depicted in Fig. 2. 
Three shown evolutions differ only with the initial eccentricity of the comet. We used 
the values ej„ = 0.01, 0.001 and 0.0001. The initial values of other orbital elements are 
ttin = 10 000 AU, ajj„ = 45°, J7„j = and ii„ = 90°. 

In Fig. 3 is compared solution of equation of motion Eq. (1) with solution of system 
of differential equations given by Eqs. (13)-(17). In both numerical solutions are the 
new term included. Evolution depicted by a solid line is obtained from solution of Eqs. 
(1) and evolution depicted by a dashed line is obtained from solution of Eqs. (13)-(17). 
Initial values of the orbital elements are ai„ = 10 000 AU, Cm — 0.3, = 60°, J?^ 
= 45° and ij„ = 45°. The Sun is located at distance 8 kpc from the galactic center, 
■^o(O) = 30 pc and .^o(O) = 7.3 km , at the time t = 0. Initial true anomaly for the 
evolution obtained from numerical solution of Eqs. (1) is fin = 0. Orbital evolutions 
obtained from numerical solutions are in good accordance. 

We used semi-major axis aj„ = 50000 AU for a comparison of numerical solutions 
of Eq. (1) and Eqs. (13)-(17) with the new terms at greater semi-major axes. Initial 
values of other cometary orbital elements and solar initial position and velocity are the 
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Fig. 6 Two orbital evolutions of a comet with initial semi-major axis 50000 AU obtained 
from numerical solution of system of differential equations given by Eqs. (13)-(17) with the 
new terms (solid line) and without new terms (dashed line). 



same as for the evolutions depicted in Fig. 3. Results are depicted in Fig. 4. Difference 
between two evolutions is caused by the large value of the initial semi-major axis of 
the comet. The initial semi-major axis is so large that the orbital period of the comet 
(T w 1.1 X lo'^ years) is comparable with the period of oscillations of the Sun around 
the galactic equator {2ti/ijJz ~ 7.3 x lo'' years). Secular orbital evolution given by Eqs. 
(13)-(17) cannot be used in this case. 

We compared also influence of the new terms on secular evolution of orbital ele- 
ments at larger semi-major axis for the initial conditions corresponding to those used 
in Fig. 4. The resulting evolutions are presented in Fig. 5. Evolution depicted by a 
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solid line is for Eqs. (13)-(17) with the new terms and evolution depicted by a dashed 
line is for Eqs. (13)-(17) without the new terms. Evolution depicted by the solid line in 
Fig. 5 is identical to the evolution depicted by the dashed line in Fig. 4. Comparison 
of Figs. 5 and 4 shows that the influence of the new terms in Eqs. (13)-(17) on orbital 
evolution is less significant than motion of the Sun which was neglected in derivation 
of Eqs. (13)-(17). 

Fig. 6 depicts numerical integration of Eqs. (13)- (17) for the case when inclusion of 
the new terms plays a relevant role. The new terms significantly changed the cometary 
orbital evolution in comparison with the orbital evolution without the new terms. Both 
evolutions depicted in Fig. 6 had equal initial conditions. Initial values of the comet's 
orbital elements were aj„ = 50 000 AU, ej„ = 0.8, Wj„ = 150°, = and ij„ = 90°. 
The Sun was located at distance 8 kpc from the galactic center, Zq{0) = 30 pc and 
Zq(0) = 7.3 km s~^, at the time t = 0. Influence of the new terms can be even more 
significant for larger semi- major axis of the comet. 

5 Conclusion 

The paper treats the effect of the galactic tide on motion of a comet with respect to 
the Sun. It turns out that the important effect from the galactic tide is the action 
of the normal «— component. The x— and y— components of the acceleration comes 
not only from the x— and y— positional components of the comet, but also from the 
0— component of the position. This is generated by the galactic disk. The effect of 
the three positional components in the x— and y— acceleration components are of 
comparable importance. 

The inclusion of the new terms into the equation of motion of the comet leads to 
orbital evolution which may significantly differ from the conventional result. This is 
true mainly for the comets with large semi-major axes. The conventional result (see, 
e. g. Fouchard et al. 2008) is obtained from Eqs. (13)-(17) putting Ti = 0, r2 = and 
^' = 0. 

The solution of Eqs. (13)-(17) is in a good agreement with the solution of the 
equation of motion represented by Eqs. (1), ifr-c(l / Tz + 1 / Tq )^^. 

We found that a comet with the argument of perihelion lu has the same orbital 
evolution as a comet with the argument of perihelion co + tt, if values of other orbital 
elements are equal. Similarly, a comet with the argument of perihelion ui and the 
longitude of the ascending node il has exactly opposite time derivatives of the orbital 
elements as a comet with the argument of perihelion tt — o) and the longitude of 
ascending node -k — O, ii values of other orbital elements are equal. 
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